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ABSTRACT 


DetmeiodcesrOrspesslole rover Movlon are defined: Synehronous preces— 
sion (without rub), partial rub, full annular synchronous rub and reverse 
whirl. A model is developed to investigate the relationships between a 
rotational system's characteristic parameters and the conditions of equili- 
brium for the full annular rub motion. The dynamics of an unconstrained 
rotor with unbalance are first established. The impact of constraining the 
rotor's motion by a rigid casing with a finite clearance is then analyzed. 
For any given system, the speed regimes in which full annular synchronous 
rub is in equilibrium are defined as a function of the system's unbalance, 
elearance, rotor stiffness, damping and the coefficient of friction between 
the rotor and casing. The rigid stator restriction is subsequently allevi- 
ated and the equilibrium conditions for full annular synchronous rub are 
reanalyzed. Additionally, reverse whirl as a possible stable mode of rotor 
motion is studied. Specific rotational systems are then analyzed according 
to the model developed to demonstrate the effect of varying particular 
Parameters. 
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INTRODUCTION 


All rotating machinery have some practical degree of unbalance which 
Gienves flexural vibrations in the system. These vibrations may be readily 
calculated in terms of the system's mass, unbalance, rotational speed, 
radial stiffness and damping imposed by the medium in which the rotating 
element operates. Should the magnitude of the vibrations exceed the 
physical constraints of the enclosure which houses the rotating element, the 
potentially destructive phenomena described here as rotor rub occurs. 

Shaft vibrations, misalignment, blade creep and thermal expansion are 
meme Ol the physical realities which precipitate design of a radial 
clearance between the high speed rotor and surrounding casing in 
turbomachinery. Theoretical efficiencies of the turbomachinery are reduced 
mieuen, by blade tip losses which result from the presence of the radial 
elearance between the rotor blade tip and the casing. Demand for improved 
system efficiency thus drives a reduction of the blade tip clearance and so 
miereases the likelihood of rotor rub. Rotor rub in fact may be a design 
condition and so accounted in the design by the use of abradable materials. 
et 

Prior to the occurrence of rub, the rotating element's dynamics may be 
described as synchronous precession, driven by the system's residual 
unbalance. If rub develops, it may occur over a fraction of the rotor's 
precessional orbit (termed partial rub), or may persist over the full orbit 
(full annular rub). Of these two conditions, full annular rub is 
potentially more destructive and is the subject of this thesis. For further 
information regarding partial rub, attention is drawn to reference [2-4]. 


Examples of full annular rotor rub occur with rotor blade-to-casing or 
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rotor-seal rub, as well as with a rotating shaft having clearance in a 
bearing. [5] The rub introduces a reaction force at the rotor's periphery 
which strongly influences the rotor's dynamics and may additionally excite 
the casing (further referred to as the stator) into motion. The reaction 
force may be broken down into its normal and tangential (friction) force 
components; the two being linearly related by the coefficient of friction. 
Friction may be moderated by proper lubrication (such as in the case of 
bearing rub) aS well as through suitable selection of materials. Other 
design considerations or abnormal operating conditions may preclude these 
preventive measures though, and so friction may take on a particularly 
destructive nature. 

Friction is the resistive force which actS tangentially at the 
interface between two surfaces in relative motion. By acting tangentially 
on the rotor's periphery, it essentially establishes a torque which counters 
the driving (input) torque for the rotating element. The result is either a 
loss in rotational speed or a nigher input torque required to maintain a 
menauant Speed. The effect of friction on the mating surfaces is to cause a 
local heat buildup and, in the very least, some degree of grinding wear and 
deformation of the rotor and seal or bearing annulus. The more extreme case 
encompasses catastrophic destruction of the machine. [2] Typically 
accompanying these mechanisms is the annoying and often unacceptable 
acouStic phenomenon due to frictional vibration instabilities known as stick 
Slip. Dynamically, the occurrence of rotor rub may destabilize the 
Synchronous precession of the rotor. If friction forces are of large enough 
magnitude, the counter-torque developed may instigate a reverse whirl 
condition wherein the precession is in a direction opposite to the direction 


meerectation of the shaft: The nature of reverse whirl is such that the 
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rotating element is subjected to high frequency alternating stresses which 


may lead to rapid fatigue failure of the shaft. [6] 
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1. DEVELOPMENT OF THE MODEL 


The characteristics of rotating machinery designed for the generation 
or transmission of power vary greatly with their particular application. An 
exact analysis of the motion is extremely complex, but many meaningful 
results can be determined through employment of reasonable simplifying 
assumptions in a single degree of freedom model. Even with these simplifi- 
cations though, numerous variables are necessary to describe the motion 
mmyolved. 

The model in this thesis is displayed in Figure 1. The rotating 
element consists of a central disk of mass, Mn mounted symmetrically ona 
laterally flexible but torsionally rigid, massless shaft of stiffness Kae 
The shaft stiffness is determinable knowing the shaft's material type, 
Meess-sectional properties and length. The shaft is then supported at each 
end by rigid bearings without clearance. The natural frequency of the 
meovor, On? the frequency of excitation at which resonance occurs, is given 
by On = 7K /M,,. The rotor disk is enclosed in its casing, separated 
radially by a symmetric clearance, Bia existing when the rotor is at rest. 
(The effects of gravity are neglected in this analysis). The casing, or 
stator, is characterized by its mass, Mo and its mounting stiffness, Ka. 
AS in the case of the rotor, the natural frequency of the stator is given by 
— 7K /M.- Damping is introduced for motion of the rotor and stator 
respectively by the damping coefficients, Dn and D. It 18 the presence of 
damping in both the rotor and stator which limits the magnitude of their 
responses at their respective resonant frequencies. Damping also causes a 
phase lag between the system's excitation and response. For ease of 


analysis, the model has been "compressed" in Figure 2 to a planar diagram 





: eee =e 
with the element's dynamic characteristics represented as in a simple linear 
resonator system: a spring-mounted mass with a dashpot in which the energy 
of the motion is dissipated. Both radial and longitudinal symmetry are 


miontecit in this model. 


-12?- 





FIGURE £ ROTOR RUB MODEL 
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2. FREE ROTOR DYNAMICS 


It is assumed that under design conditions, the rotor is driven at a 
constant speed, w, independent of the external load applied. This assump- 
tion essentially neutralizes the effects of fluctuations in the design load 
due to damping or rub related friction. Under ideal conditions, it is 
imaginable that the rotor is perfectly balanced. Thus when driven at its 
design speed, the rotor center remains at the center of rotation; aligned 
with the axis defined by the supporting bearings’ centers. This idealiza- 
tion is not realized in practice though, and to some extent, a residual 
unbalance exists in the mass distribution of the rotor. This unbalance may 
be depicted in the planar model by locating the mass center (or center of 
gravity), g, at a distance 6 from the center of the disk. The effect of the 
unbalance is to establish a centrifugal force acting upon the center of 
gravity while the rotor is in motion. This centrifugal force, which is 
proportional to the square of the rotational speed, acts as the driving 
moree Tor the rotor's dynamical motion. 

An analysis of the motion of the free rotor (i.e., without displacement 
constraints imposed by the stator) provides the baseline condition for which 
the rub condition may be compared. The equation of motion for the free 
rotor is derived by applying Newton's 2nd Law to the rotor: 

(sum of forces acting on the rotor) = 

(mass of the rotor) x (acceleration of the rotor center of gravity). 
Referring to the free body diagram for the rotor, Figure 3, the forces which 
must balance in the steady motion of the rotor are described in the 


following: 
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FIGURE 3. FREE-ROTOR FREE BODY DIAGRAM 
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a. Centrifugal Force - As described above, this force results from the 
presence of a residual unbalance; i.e., the center of gravity does not 
coincide with the center of rotation. The centrifugal force acts on the 
center of gravity and is directed radially outward from the bearing axis 
(stationary reference position about which the rotor precesses). 

b. Spring Force - The spring (or restoring) force acting radially 
inward on the rotor center results when the rotor center is displaced from 
the bearing axis. It is assumed in this model that the Spring force iS a 
linear function of the shaft stiffness (assumed constant) and magnitude of 
displacement. 

ea Daiiping Force — The damping force acts on the rotor center such 
that it resists translational motion by the rotor. The damping force is 
assumed to be a function of the damping coefficient experienced by the rotor 
in its particular medium and the linear component of the rotor's velocity. 

The nature of the rotor's dynamics readily lend the analysis to 
description using complex variables. Defining a as the time dependent 
Bosition of the rotor's center, referenced such that Z. = 0 coincides with 
the rest condition wnere the rotor center is aligned with the bearing axis, 


then the equation of motion may be written as: 


Po Meee eee 
- mae rée7 neg (1) 


> 
where Z ait) = the vector position of the center of gravity of the rotor. 


> 
Furthermore, Ae may be broken down into the vector sum: 


> 
where age Geplede MUSmene VCCLOr Position of the rotor center of gravity 
melative to the disk center. 


The equation of motion then becomes: 





M ( ~ D > K as 
+ = - — 
iG ae 2 30 rc | Gea 
or, 
a a > > 
Maze baz KZ = fz (23) 
rc re re rec 


Zoo (t) is simply a function of the radius of unbalance and rotational 
speed: 


z co sala 
ge 


Assuming synchronous whirl response, the motion of the rotor may be 
described by: 
> > jwt 
z(t) = |zle 


Making these substitutions into Equation (2), and simplifying, results in 


the following: 


This harmonic motion of the rotor is termed forward synchronous precession 
peice the orbital motion of the rotor is characterized by the forcing 
function's angular frequency and direction. Breaking down the rotor's 


mechanical impedance, 
[(K, - wM,) + j(uD_)] 
r r ic ; 


into its magnitude and phase then leads to the harmonic response equation: 


= 
Zz, iG a ee ee we 
a ) C3} 


The phase angle ¢ describes the lag between the rotor's response and the 


forcing function. This phase lag which is due to damping is given by: 
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It is convenient to introduce non-dimensional variables into Equation 
(3) for ease of future analysis. By defining a non-dimensional rotor speed 


and damping ratio as given below: 


> 68 J (wt->) (4) 
acieenety = (ca) 


Similarly, the phase relationship may be rewritten: 


¢= tan -*222+=- (4a) 


The plots of Z vs. B and 9 vs. 8 are typical amplitude and phase plots 
for a second order system, such as those in Figure 4. 
Consider tne case of light damping, c << 1.0. At low Speeds, i.e., 
B << 1.0, the magnitude of Equation (4) reduces to iz. = eae As the speed 
increases, the displacement increases non-linearly until the maximum dis- 


5 Wee 
placement occurs at £ = (1 - 2z7°) 


— a ee ee ee ee ee ee oe 


Here it can be seen that the maximum displacement is solely a function of 
the radius of unbalance (directly) and the damping ratio (inversely). 
Further increases in speed bring a reduction in displacement. For very high 


speeds, i.e., B >> 1.0, the displacement approaches 6 (a constant). 
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Tne phase relationship for light damping is given by Equation (Ha). 
For low speeds, 8 << 1.0, 6 iS approximately zero. Essentially, the forcing 
function and the response are exactly in phase. As the speed increases 
though, and approaches the resonant speed, the damping force becomes 
Significant and the response begins to lag the forcing function. Near 
resonance, a rapid phase shift occurs such that at resonance, » = 90°. As 
the speed further increases, $9 asymptotically approaches tm. In the limit, 
ie., 10r © = 0; 6 = O for all values of 8 < 1.0 and $ = wm for all values 


orp 7 1.0. 
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3. CONSTRAINED ROTOR DYNAMICS 


As previously described, rotating machinery designed for the generation 
or transfer of power are typically enclosed within a casing which minimizes 
the clearance between the rotating element and casing. The dynamics of the 
constrained rotor motion may simply reduce to those of the free rotor as 
Biesented in the preceding section if the rotor and stator never come into 
Memeact. This no-rub condition exists if either of two basic criteria are 
ma effect: 


(1) The maximum amplitude of free rotor displacement, is less 


|2aleees 
c'max 

than the constrained rotor clearance, or 
(2) The rotor is operated within a speed regime, characterized by a 


> 
Speed w, such that the amplitude of free rotor displacement, jz, (w) |, is 


less than the constrained rotor clearance, 


As will become evident through the following analysis, condition (2) above 
should be further restricted to speed regimes where w is less than the rotor 
natural frequency, Won Tnese two cases are depicted in Figure 5 uSing the 
plot of amplitude vs. speed. 

If in fact the rotor and stator come into contact, the dynamics may 
become quite complex. The analysis will be first conducted with the 
assumption that the stator is perfectly rigid (though this restriction will 
be alleviated later) and only full annular rub will be considered. 

The free body diagram for the rub case is presented in Figure 6. Prior 
to contact, Equation (4) describes the rotor's motion. It can be Seen that 
Meoen the initiation of rub, the rigid stator limits the rotor's displacement 


° . a ° a * ° 
While introducing a reaction force, R(t), at the point of contact. This 





ROTOR DISPLACEMENT 


ROTOR DISPLACEMENT 


(Zc [vc) 


(Zc/tc) 


e 


OF I 5 Se aoe 
WO-NUAUDNAGO 


9909 
aN O 


0.5 
0.4 
0.5 
GiZ2 
Ox 

0.0 


O —-* —- —* 8 -8 -*# 8 -* -* os fj 
OoOoO-NUbAYUOnN @WOO 


FIGURE SCa). 
FIGURE SCb). 


== 


Rotor Displocement vs Speed (No Rub) 
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FIGURE 6 ROTGR RUB FREE BODY DIAGRAM 
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reaction force is comprised of a normal force and a tangential (friction) 
force. Since the displacement is constant, the restoring force acting on 
the rotor is also constant (i.e., independent of rotor speed, w). As in the 
ease of the free rotor, the equation of motion for the rub condition may be 


developed from Newton's Second Law: 


> > 
z= De Khe ee en (6) 
ome re re 


, > 2 > > jut 
Substituting oe = (Zz + Z)s and R(t) = (N + juNj)e , then Equation 


(6) becomes (after rearranging terms): 


5 > a 
Mz+Dgzg+Kz =-Mz rt) 
re re re r° ge 


Since Z oc = ce oe then the above equation may be solved for Z. by 


assuming synchronous precession: 


[Mw - (N + §uN)] 
ZA) eee Se e (7) 


By dividing through Equation (7) with the radial clearance, ry» and 
utilizing the previously defined non-dimensional speed and damping ratio 


variables, Equation (7) is transformed to its non-dimensional form: 


Eom ep (oun) ye aJut (8) 
: [(1 - g°) + 5(278)] 


where the new variables are defined: 


e = 6/r 
C 


n N/kr 
C 


The right hand side of Equation (8) is reduced to a magnitude and a phase in 
the following rotor rub equation: 


1/2 
ieee n)° “ Ghaal JC wt 2 ny SO) 


- > 2S ie (9) 
¢ uC) = (Gyo eee nies 





= 5 = 


where 


The following observations may be drawn regarding the physical 
iierr1cations of the rotor rub equation. The parameter 6 1S a measure of the 
rotor's unbalance, and as such is not designed into the rotor, but results 
from manufacturing processes which are unable to meet with ideal design 
conditions. The value of ce though is designed into the system, and so 
tolerances muSt be placed on 6 in order to ultimately control the value of 
€, a measure of the rotor's displacement-to-clearance ratio. In the case of 
m@emiree rotor, any increase in 6 leads to a direct increase in the rotor's 
maximum diSplacement. Since the maximum displacement is now constrained 
though, a normal force develops at the interface between the rotor and 
stator when the rotor's runout equals the clearance. The maximum normal 
force developed thereafter is directly related to the value of e€. Referring 
memrieure 5, the normal force developed for rotor speeds less than the rotor 
resonant speed may be crudely estimated as the magnitude of force necessary 
to compress a ‘Spring' of stiffness equivalent to the rotor's impedance a 
distance equal to tne difference between the free rotor's Concent and 
Maeeradial clearance. For speeds greater than the rotor resonant speed 
(8 > 1.0) the characteristics of the combined rotor-stator govern the 
system's behavior and a more detailed analysis is necessary. 

Meanwhile, a frictional force, f = wn, develops at the rotor-stator 
interface in opposition to the relative motion between the rotor and stator. 
Meescocfficient of friction, u, depends primarily on the material types in 
Bee@edet and the degree of lubrication present. For unlubricated metal-to- 


metal applications, typical values of yw may range from 0.1 - 0.3. [7] The 
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frictional force may be replaced by a force and a couple acting upon the 
rotor's center. The torque thus created acts as an additional load upon the 
rotor's driving force. Equation (9) emphasizes the phase lag, ($+ Y), 
between the response of the rotor and the forcing function. 9 was defined 
previously as the phase lag caused by the damping force in opposition to the 
linear component of the rotor's velocity and was given by the expression, 
> = tan | (228) 70 = ey Y is the phase lag caused by the frictional force 
which further opposes the rotor motion as previously given by: 
Y= tan | (un)/(e8* = Aig 

Since rub iS presumed, the magnitude of the rotor's displacement is 
exactly equal to the radial clearance (given a rigid stator). This 
condition results in the magnitude of Equation (9) being equal to unity. 
The independent variable in the equation is the rotor speed. The values of 
u, © and e€ are assumed constant, and so the values of the normal force and 


phase angle may be solved for any value of g. This arises from the 


mend ition: 


Meraen Giny( |< 


exe a 
Meee bo) ete Cece) | 
Following a few algebraic steps, the following quadratic equation inn is 


derived: 


2 2 ene 2 4 
(1+ uo)n© + 2(1 - Bo + Qcgu)n + [1 - Boye + (208)° - eB 1 = 0 ~~ (10) 
Which when solved yields the relationship for nas a function of 8: 


|e 


ee ee = Se Se 6 eS Se ee ee = Se Se Se Se ee Se Se eee ew tee we een ee SS 


ye 


Equation (11) indicates two possible solutions for the normal force. 
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fmemisebie Upper, ‘positive’ root which is of primary interest, yielding a 
condition of n = 0 which is consistent with the rotor speed, 8, at which rub 
initiates. As 8 varies, so too the values of the normal force and phase 
angle vary together as necessary to maintain a balance of all forces acting 
eaeche rotor. 

It is evident from Equation (11) that throughout the regime of possible 
rotor speeds, depending upon the values of the system parameters u, ¢ and e€; 
the normal force may be positive, negative or imaginary. A negative value 
of the normal force here implies tensile stresses existing between the rotor 
and stator. This is essentially infeasible, and so for all values of n < 0, 
the rotor and stator are deemed not to be in contact. The speed at which 
rub is actually initiated is determinable by solving Equation (11) for the 


Value of 8 at which n = 0: 


Mais value of 8 then results in Iz. equal to r, in both Equations (4) and 
(9). It should be immediately noted that if e/2c < 1.0 for a particular 
Seeevem, then rub in fact does not occur at any rotor speed. The inequality 
above is simply the numerical equivalent to a previously described condition 
- i.e., the maximum amplitude of the free rotor displacement is less than 
the radial clearance allowed the rotor. This is indeed the most common 
designed-for condition. Through service use though, the value of e for a 
rotating machine may expectedly change since both the degree of unbalance 
and the radial clearance are affected by such physical realities as creep 
caused by high centrifugal loading, thermal and aerodynamic stresses caused 
by interaction between the rotor and the working fluid, or possible shaft 


bowing due to uneven heating or cooling or simple static loading. 
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It is the intent of this analysis to determine the behavior of rotating 
machinery as the rotor speed is increased beyond the rub speed for systems 
where ¢/2c > 1.0. Figure 5(b) provides a simple plot of the amplitude vs. 
Speed relationship for the rotor. It is known that for B < eee Equation 
(4) describes the rotor's motion. Furthermore, for rotor speeds within the 


regime Jae < 68 < B', it is assured that the rotor and stator will be in 


Db 
contact with the normal force given by Equation (11). For rotor speeds 
greater than g' though, Figure 5(b) can only provide basic insight into the 
problem. It was alluded earlier that while the rotor is undergoing a full 
annular synchronous rub, then it is the characteristics of the combined 
rotor-stator which govern the system's performance. This performance may 
best be analyzed by investigating the normal force and phase angle 
experienced by the rotor. 


Equation (11) for the normal force may be rearranged algebraically to 


yield the following expression: 


Wee 





2 nes goes ee oe eo SO See ee en en Bae (12) 


The three possible regimes of n (positive, negative or imaginary) may be 
examined through this equation. it is instrumental to first analyze F, the 
Sum Oi the external forces acting on the rotor center, by resolving it into 
Orthogonal components (see Figure 7) - those forces acting parallel to om 
(designated the radial or r-direction) and those forces acting normal to 
ze (designated the tangential or t-direction). 

F is further resolved into a component due to rotor displacement and 
eeeomponent due to the mass eccentricity itself in Figure 8. 


By combining the non-dimensional relationships in Figures (7) and (8), 
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pmeeonetric interpretation for the relationship between forces in Equation 
ay) iS attained in Figure 9. 
By applying the Pythagorean Theorem to Figure 9, the variables become 


related by: 


(e8°)* = (2¢8 + un)© + [(1 - Bo) + nn]? (13) 
Mayeiecally, this states that: 
meee component of centrifugal force acting on the rotor center due to mass 
eccentricity], ee, is balanced tangentially by [the damping and friction 
forces], and radially by [the normal and elastic forces lessened by the 
centrifugal force acting on the rotor center due to rotor displacement]. By 


solving for the normal force itself, Equation (13) yields: 


yee 
eee ees) > — (ere + uin)*} (14) 


The parallel between this equation and Equation (12) is immediately 
apparent. The following generalizations may now be drawn regarding the 
nature of the normal force as a function of the remaining system parameters 
Emaevariables: 

(a) As the rotor speed increases: 


(i) The normal force, n, first becomes positive when: 


Wee ec) (ere) 


This simply states what was earlier determined; that rub initiates when the 
magnitude of the centrifugal force due to mass eccentricity equals the 
rotor's impedance multiplied by the radial clearance. 

(ii) The value of (8° - 1) increases; which tends to increase the 
normal force. 

(iii) The value of the centrifugal force, Ce increases; which will 
increase the normal force if Uae phase angle, 6, is less than 90°; and will 


decrease the normal force if » exceeds 90°. 
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(iv) The value of the damping force, (278) increases; which drives the 
phase angle towards 90°. For values of » less than 90°, this means a 
reduction in normal force, and obversely, for values of » > 90° this 


translates to an increase in normal force. 


(b) Any increase in the damping ratio, zg, or friction coefficient, u, will 
increase the tangential forces, and in a manner similar to (iv) above will 
drive the phase angle towards 90°. Again, for values of ¢@ less than 90°, 
this means an increase in normal force, and for values of @$ greater than 


90°, this means a reduction in normal force. 


(c) The center of gravity of the rotor will orient itself to that phase 
angle necessary to balance the forces acting upon the rotor. This assumes 
that there is in fact some angle o which will satisfy this dynamic 


requirement. 


If for a given rotating system characterized by the parameters u, zt and 
€ operating at a particular speed £; there is no value of the phaSe angle 
which will cause the forces acting on the rotor to balance, then the 
assumption that the rotor is precessing at synchronous speed is no longer 
valid. When this case arises, i.e., when the radical term in Equation (12) 
becomes negative and n becomes imaginary, an alternative mechanism for the 
rotor motion must be determined. The ‘equilibrium criterion' for 
synchronous rub may thus be stated mathematically: 


Iam 


me 2 
(ep°)° = DiGi Po Cece) i a 0 (15) 


then synchronous rub is not a possible mode of rotor motion. 
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The multiplicity of rotor speeds (roots) which make the radical term in 
Equation (15) identically equal to zero lead to three possible Sequences of 


rotor modes of motion occurring within the rub regime: 


(i) If all the roots to Equation (15) occur at speeds below the speed 
at which rub initiates, then a real solution to the rub equation exists for 
all rotor speeds. This does not preclude the possibility that rub may cease 
(i.e., n< 0), but rather indicates that the rotor dynamics will be governed 
by the relationships previously developed for synchronous whirl (i.e., 
Equations (4) and (12)). 

(ii) If one of the roots to Equation (15) occurs at a speed greater than 
the rub speed, then the rub equation has no solution for rotor speeds 
greater than this limiting speed Br: It is possible that BT Could eccur at 
a speed for which rub is not in effect (i.e., the normal force is less than 
zero according to Equation (12)). This condition shall be discussed in 
greater detail in a later section. 

(iii) If two of the roots to Equation (15) occur at a speed greater than 
the rub speed, then there is a limited regime within which synchronous rub 
is not in equilibrium. For all rotor speeds less than the first root and 
greater than the second root, the rotor dynamics may be described by 
Equations (4) or (12). Otherwise synchronous rub is not a possible mode of 
meror motion. 

Throughout the speed regime in which Synchronous rub is not a solution, 
three modes of rotor motion are possible - either a forward synchronous 
whirl without rub, a partial rub or a full annular rub reverse whirl. These 
cases shall be discussed in greater detail in the following section. In 


order to provide clarity on the interrelationships of these various modes, 
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emewoutline in Figure 10 provides a sequencing of the possible system 
Miiamics) as a function of rotor speed. It is informative to demonstrate 
these various modes by comparing the amplitude of displacement and normal 
force experienced by the rotor as functions of rotor speed as in Figures 11 


through 14. 


oi) No rub 

The reference case is the simple condition that rub never occurs. This 
is true for systems where e« < 2t. AS can be seen in Figure 11, the radial 
clearance is greater than the maximum displacement experienced by the rotor 


and so the normal force is zero for all speeds. 


eo.) Limited synchronous rub 

As depicted in Figure 12, rotor rub initiates at some speed ae oie 0. 
Further increases in rotor speed drive the rotor hard against the stator 
wnile the normal force steadily increases. The large friction and damping 
forces which also develop ultimately disrupt the synchronous rub. Further 


me@emeases in rotor speed do not lead to a restoration of rotor rub but 


rather follow the course of action for the free rotor. 


ee) Synchronous rub-limited discontinuity 
Similar to the limited synchronous rub case, rub initiates at a speed 
8 


< 1.0, and ceases at some speed, 8 hee By» rotor precession jumps 


nub lt. 
down to the free-rotor synchronous whirl. As shown in Figure 13, at a yet 
higher speed, dynamic conditions once again become favorable for a stable 
Synchronous rub. At this speed though, some perturbance is required to 


bring the rotor into contact with the stator and thus initiate the rub. It 
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must be emphasized that reverse whirl may also be stable at this Speed 
(dependent upon the slip condition), thus making it difficult to predict the 


exact motion which would ensue. 


(iv) Continuous synchronous rub 

Pemtmerne limited rub case, Trovor rub initiates at some rotor Speed, 
Sb < 1.0. As rotor speed is further increased, the rub reaction force 
also steadily increases. In this case though, the damping and friction 


forces are not significant enough to disrupt the equilibrium of the 


synchronous rub (see Figure 14). 


pmalyiically, regimes where the varlous modes dictate the behavior of 
the rotor may be defined and mapped onto a family of curves using variables 
malioeiar defined. Figure 15 depicts the impact on rotor motion of varying 
the coefficient of friction, unbalance (or inversely the radial clearance), 
Or speed ratio for a fixed value of rotor damping. For this particular case 
@eeee= .10, region (i) of the plot depicts the condition for which no rub 
occurs (independent of the rotor speed) and is described by the inequality 
Meee < 2c. For a particular value of « > 2z, rub first initiates at the 


eweea defined by the previously derived relationship: 


2 iy a am at = a 


This value provides the lower bound to region (ii) in Figure 15. Region 
(ii) depicts the speed regime wherein synchronous rub is in equilibrium. If 
we consider the specific case where the coefficient of friction, yp = 0.25, 


then region (iii) depicts the speed regime where synchronous rub has fallen 
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oeeor Caquilibrium. The lower bound to this region 1s defined by the 


inequality provided in Equation (15): 


2.2 
cco toa. = SS a =" 0 
Clr aie) 


Synchronous rub never re-establishes equilibrium in rotational systems with 


a value of e described by the following inequality: 


Otherwise, the speed at which synchronous rub recurs is depicted by the 
curve labelled "yp = 0.25." This curve is generated similarly to Equation 


(15) above, but now represents the case: 


(a a 
cegrye - bul r_B)_+ (2zB))° . 5 


Cie 1 


Lastly, for any system described by g¢ = 0.10 and yp = 0.25; if 


1/2 
e 2 (2° ee ae. + i) i, then synehronous rub is a Solution to the rotor 


motion for all speeds, 8 > Bub’ A plot such as in Figure 15 thus clearly 
Outlines rotor motion as a function of speed for any system described by the 
variables uw, c and e. Specifically, each of the four cases described in 
this section are delimited according to the system's value of e: 


Mee no rub: « < 27 


(ii) Limited synchronous rub: 


mer) Synchronous rub - limited discontinuity: 


- & x uo 
ee 


ea) ale ia) 





zy Pe 


(iv) Continuous synchronous rub: 


The relationships above explicitly indicate the effect on the rotor's mode 


of motion of altering any of the variables involved. 
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4, DYNAMICS OF REVERSE WHIRL 


Reverse whirl is the state of rotor dynamics wherein the direction of 
precession of the rotor center about the bearing axis opposes the direction 
of rotation of the rotor about its own center. It was shown that in the 
absence of rub, an unbalance in the rotor's mass distribution results ina 
forward synchronous whirl of the rotor. The presence of damping resists 
ears tOrward precession and thus creates a phase lag between the rotor 
center as it precesses and the mass center as it rotates. This damping 
foree though, is linearly related to the whirl speed and so can not by 
itself destabilize the precessional motion of the rotor. When rub occurs, a 
frictional force at the interface between the rotor and stator develops 
wnich combines with the damping force to oppose the rotor's forward 
M@eeession. Unlike the damping force though, the frictional force is a non- 
linear function of rotor speed. If the additional force created by the 
frictional force opposing the rotor's forward whirl can not be balanced by 
the tangential component of the centrifugal force associated with the 
rotor's mass unbalance, then the synchronous whirl assumption breaks down. 

The dynamics of the reverse whirl problem are far more complex than for 
forward synchronous whirl. Since the whirl speed, 2, is now independent of 
the rotational speed, w, it must be introduced as a new variable in the 
analysis. This then creates the additional complexity of having a phase 
angle which is now a function of time (more specifically, a function of 
(w - 2)t). Since the phase angle varies regularly with time though (i.e., 
constant frequency), a simplifying assumption is that the time averaged, or 
net position of the mass center of the rotor is located at the centroid of 


the rotor (i.e., 6 = 0). Furthermore, since the whirl and rotational 


ee 


directions are opposed, then the direction of relative motion at the point 
of contact (designated point A) between rotor and stator must be determined. 


The velocity of point A may be written: 


<= 
WW 

<1 
i. 


= 7 rf = linear velocity of the rotor center 


|v = aw = linear velocity of point A relative to the rotor center 


Ao 
where a = rotor diameter 

By assuming (aw) > (r2), then the linear velocity of the rotor at A may be 

assumed in the direction of the rotor's rotation and the friction force is 

thus aligned with the direction of whirl. 

The free body diagram may now be analyzed to determine the conditions 
whereby reverse whirl is possible. Figure 16 displays the forces at work 
and their assumed orientation for reverse whirl. By considering the balance 
Of forces in the tangential and radial directions, two independent equations 
are developed. In the tangential direction, the friction force is balanced 


solely by the damping force: 


uN D (ra) 
By defining p = 2/w and substituting non-dimensional variables: 


eTBp (16) 


un 
In the radial direction, the centrifugal force of the rotor is balanced by 
the elastic force of the rotor plus the normal force between the rotor and 


meacor: 
Z 
Mr -=N+Kr 
an ve. Ee 
In non-dimensional form this becomes: 


note eo (17) 


The rotor rotational speed, w, and whirl speed, 2, may be kinematically 
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=15— 


related through the slip condition that exists at the rotor~stator 
interface. A no-~slip condition would yield the kinematic relationship 
common to epicyclic gear trains: 

Q = (a/r dw 
If there is slip at the interface, the whirl speed is reduced accordingly. 
The slip ratio, s, is thus defined as the ratio of the (rotor whirl speed)- 


fe-the- (rotor whirl speed without slip); i.e., 


Equations (16) and (17) may thus be rewritten: 
ie 2cbl (a/r)s] (16a) 
[s(a/r,)1°B° = n+ 1 (17a) 
By substituting (17a) into (16a), the critical condition for reverse 
marl is determined: 
wiL(a/r,)°s°1B° - 1} = 2¢6l(a/r)s] 
The above condition may be resolved into a quadratic equation in 8: 
ul(a/r)s]°B° - 2cl(a/r.)s]8 - u = 0 (18) 
Solving for 8 yields the critical speed above which reverse whirl is 
possible: 
ee 


a 2c 
ane eect) 
ee Ge alt aii wa oo J (19) 


The lone unknown variable in the above equation is the slip ratio, s, 
wnich may in fact be difficult to determine. The limiting case arises of 
course when no slip is present (i.e., s = 1.0). If it is assumed that the 
damping and friction are of the same order of magnitude, then Re is of the. 
order of ri/as wnich in general iS much less than unity. The limiting case 


thus states that reverse whirl is possible for esSentially all rotor speeds 
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mieinterest (1l.c., Witmin the rub regime). It is not until a significant 
degree of slip is present (i.e., s is of the order of r fa) that the minimum 
possible speed for reverse whirl may become restrictive on the behavior of a 
Recor undergoing full annular rub. 

Since the rotor first assumes a forward synchronous precession though, 
if will continue with this motion until either it becomes unStable or is 
perturbed into another stable mode of motion. As previously developed, the 
morward Synchronous whirl without rub 1s continuously stable as long as the 
displacement of the free rotor is less than the constrained rotor's radial 
clearance. The Only mode where synchronous whirl may be destabilized then, 
is during a synchronous rub. This has been analytically developed and 
occurs if the reaction force between the rotor and stator becomes imaginary 
or if the equilibrium state for the rotor-stator is inherently unstable. 

In this situation, the equations developed based on synchronous whirl are no 
longer applicable. 

A qualitative analysis of the transient state following destabilization 
of the synchronous rub proves quite meaningful. It may first be reasonably 
presumed that the rotor will assume either a partial rub or a forward 
subsynchronous whirl following the destabilization. Considering the latter 
mode of motion, the impact of the non-synchronicity on the rotor's dynamics 
is primarily that the phase angle, 9, iS now a time dependent variable. The 
centrifugal force which has been pressing the rotor against the stator 
throughout the full annular rub now has a fluctuating component. As the 
phase angle approaches 180°, the centrifugal force acting on the rotor 
eemeer due to tne mass unbalance is actually joining with the spring force 
to drive the rotor center in towards the bearing axis. The normal force 


rapidly drops to zero and contact is broken between the rotor and stator. 
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This may be somewhat anticipated by observing the near-vertical negative 
Slope of the corresponding normal force plot at the critical equilibrium 
speed (see Figure 17). Now that rub is no longer present, there is no force 
to drive a reverse whirl, and so the rotor will jump down to the stable mode 
©f forward synchronous whirl without rub. Perturbing the rotor back into 
contact at this speed (or higher) may now provoke a reverse whirl. In the 
brevious Section, the possibility for a restabilization of synchronous rub 
was demonstrated. Under these conditions, either synchronous rub or reverse 
whirl may result when the rotor is perturbed back into contact with the 
Stator. In general, no smooth transition can occur from a synchronous whirl 
to a reverse whirl since this would mean transiting through a regime of zero 
whirl speed. The whirl speed though would not drop below that necessary to 
maintain rub-contact; for when contact is lost, the rotor either assumes a 
partial rub whirl or resumes a synchronous whirl without rub. 

The behavior of the rotor when transiting from a higher to a lower 
speed differs markedly from the behavior of the rotor as it is increasing 
speed. In the absence of some outside perturbation, when reducing speed 
from a Situation of synchronous whirl without rub at a speed, B > 1.0; no 
jump phenomenon in the rotor displacement curve would occur, but rather rub 
would occur only within that region of Figure (5b) where the free-rotor's 


displacement exceeds the rotor radial clearance. 
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5. ROTOR RUB WITH A NON-RIGID STATOR 


The assumption of an ideally-rigid stator enclosing the rotating 
element in this analysis is not feasible, but depending on the stator's 
eynamice Characteristics, such may be an entirely satisfactory assumption for 
the analysis. In this section, the stator is assigned a finite mass, radial 
stiffness and damping coefficient. For ease of analysis, these characteris- 
tics are assumed to be axisymmetric. Additionally, the stator is assumed to 
be torsionally rigid. It is often instrumental to investigate the 
translational system analogous to the rotating system under analysis. This 
is first accomplished to provide insight to the complex equations developed 
in the final analysis. 

The basie precept in analyzing the coupled rotor~-stator system is that 
once the state of continuous contact is interrupted, the synchronous full 
annular rub condition will rapidly deteriorate to a partial- or no-rub 
Seate; or if sufficiently perturbed, to a reverse whirl. The translational 
system designed to model the rotor-stator is presented in Figures 18(a) and 
(b). The rotor and stator are assigned their respective masses, a linear 
Spring and damper to characterize their dynamic characteristics, and are 
restricted to motion along the horizontal axis. In the rotational system, 
Steady motion of the combined rotor-stator may be depicted by an are which 
traces the path of the rotor-stator center of mass. The actual force 
experienced by the combined mass would have a frequency equal to the rotor's 
precessional frequency, anda magnitude which depends upon the centrifugal 
force acting on the rotor and the dynamic characteristics of both the rotor 
and stator. In the translational system a static pre-load, Po? is presumed 


to displace the rotor mass deross tne tadial elearance, oe and then further 
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displaces the combined rotor-stator mass to a static equilibrium position 
moouc which the system vibrates. A sinusoidal force then drives the system 
at the 'synchronous' speed, w (analogous to the rotor's rotational speed), 
with a magnitude Py. 

The free body diagrams for the rotor and stator under the static pre- 
load are displayed in Figure 18(c). 

The force which actS at the interface is given by Pos Essentially, 
only the stiffness of the rotor and stator resist the static load. The 


Static equation for the rotor and stator are given by Equations (19) and 


(20) respectively: 


Pw eee (19) 

ES = K 5X5 (20) 
Additionally, X, and X. aresretaved=by the kinematic condition: 

* - OG * “ 


Therefore, a single equation may be written: 
_ = ~ 
ES in k,(r, P /ks) 
Solving for the static load at the interface yields: 
a oe = Kor kA /{k, : k,) (21) 
The equation of motion for the system's dynamic condition is similarly 


determined by analyzing the forces acting on the rotor and stator. The 


dynamic condition free body diagrams are provided in Figure 18(d). 


Applying Newton's 2nd Law, the balance of forces may be written: 


mOvor : 
Be Fe . _ jwt is jJwt : ve 
KX, DX, - Pie F Fe mx, (22) 
SaLlor : 
= _ : jJwt _ Ac 
kK 5X5 Dx, € PS cf Pie = M,X, (23) 
Where P Boe Eopaeceilcomeunenaymamic reaction force at the rotor-stator 


1 





‘ 
‘ 


_—s 


=o |= 


Vi ay ia ae 


F 
xy | Xp 
D, Bb. 
| 
atl mp 
ROTOR STATOR 
N i 
kK | J Kp 
€ 
' D 
< a 2 es 
~ 
Se 
~~ 
VN ky kp 
Fo 
kK, x ry Pe Po Ke Xp 
Fo +f; eJut 
Ke Xp 
aes 
P, +P, eJut —ere De Xe 





FIGURE 18. TRANSLATIONAL RUB MODEL 
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mrceri ace. 
By aSSuming the coupled system's response is a harmonic motion of magnitude 


A, and recalling the kinematic relationship between x, and x the rotor's 


] Be 


and stator's displacement variables may be rewritten: 


jut, 


as 
il 


+ 
oar PK + Ae 


jut, 


X 


. oes + Ae 


2 


Making these substitutions into Equations (22) and (23) results in the 


following equations of motion: 


BOLor: 

a JGR : JOC jJwt Je 2 Jue 
k, (r +P o/kotAe ) D, (jw) Ae Des Pie * ae +r Fe = m, dw 2 
stator: 
- Jwt,) _ Jwt Jut _ eine 

k,(P (/k,* he ) D( jw) Ae + oe + Pie = mw e 


These relationships are valid only if the rotor and stator elements maintain 
GOntinuous contact. Eliminating the static equilibrium relationships from 
the above equations and solving each equation for the dynamic reaction at 


Mmiesincverface yields equivalent relationships: 


2 

P F + At k, jwD, + mw ) 
2 

Py 2 7 mu) 


Roe 
ACK, jwD 
These two relationships may be reduced to a single equation and solved 
for the magnitude of the vibrational response: 
BRE ans roe 5 OS ae or ean (24) 
[(k,+ k) ~ w (m+ m,)] + jw(D,+ D,) 

The criterion that 'synchronous rub' (as portrayed by this model) 
remains in equilibrium is evidenced through the model by a positive 
reaction force at the interface, i.e., that the magnitude of the dynamic 


reaction force be less than the static reaction force: 


TP] S19, 
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Scoipimine sauatvions (22), (23) and (24), this criterion becomes: 
2 : 
AL(k,- wm,) + jwD,I]| < (Fo~ kyr )k,/(kyt ky) (25) 


|F, 


2 
a emote eel ee > kt )kK5/(ko 4 Kk.) (26) 
Substituting for the variable A using Equation (24), this condition becomes: 


2 ; 
FUCK, w) m,) + jwD, J k 
1 


It is convenient to employ the following definitions when referring to the 


rotor, stator and combined rotor-stator: 
wig) = ECs Gas m)1'/* 
Or, = (k,/m,)'”? 
Wo = (k,/m,)'"* 
Sh re D,)/2L(k, + k,.)(m, + m)1!/* 
6, = D,/2ck,m, 11" 
ee D,/2L km!’ 


The condition that "full annular synchronous rub" be a solution to system 


motion now becomes: 


2 
oa ee oe 
; L(1 2 ) j(26, i a : 
ee er Le tN yn INE 8 - oe 
Py Paks, eC eo oS Roe Me 
Li = ey, + J6 25. =a) al 
Oy N 


Tnis inequality may be further simplified by eliminating common factors and 
by writing the two complex dimensionless impedances in terms of their 


magnitude and phase. The phase relationships which evolve may be 
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disregarded since only the magnitudes involved affect the "equilibrium 
eriterion". Rearranging terms results in the final criterion dictating the 


equilibrium of synchronous rub: 


[ae (F - k,r ) 
ee ag 
eon = “75) + (20, ie ] Cae = ware a (20, =e sae 

ON Bile AN 2 
The inter-relationships of the various terms: |e relative to 
2 2 2 Pale : 
(F kr); and td (w / Uy i + [2¢ (w/w ) J } relative to 
Ei Z Fa CA oP : P : 
{[1 (w / Wn 5 ) Se se [2¢,,( w/w.) J } outline various regimes in the 


fetereoi (force vs. frequency) where synchronous rub is in equilibrium. Four 
fundamental cases are depicted graphically in Figures (19) through (22): 
Sei kar) > |F,{ and uw. > wy 
eee re kr) < «(Fj and uw. > uy 
ee (Fo - kr) > |F,| and Us < Wy 
Sipe (Fo kr.) < fae and uw. < wy 


Throughout tnese figures, the shaded region indicates the speed regime 

wherein synchronous rub does not satisfy dynamic equilibrium in accordance 
with the condition provided by Equation (26a). The artificialities of the 
translational model are evidenced by the awkward low frequency behavior of 
the graphs. In particular, there is no clear correlation between the terms 


(FO ce kr) and F and similar functions in the rotational system. A 


1? 
Single most prominent observation in all four cases though, is that 
ememronous rub 1S not in equilibrium in the vicinity of the combined rotor- 
stator natural frequency. This supports earlier statements that it is the 


characteristics of the combined rotor-stator which govern the stability of 


synchronous rub. 
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FIGURE 19 RUB CONDITION - CASE 1 (Top) 
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FIGURE 22. RUB CONDITION - CASE 4 GBottom 
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The more exact analysis of the rotor-stator dynamics requires the 
assignment of finite mass, radial stiffness and damping characteristics to 
the stator in Figure 2. As earlier stated, these characteristics are 
assumed axisymmetric. The reference position for both the rotor and stator 
displacements remains the imaginary bearing axis, B, which marks the rest 
Mesrcion for the rotor and stator. Figure 23 depicts the rotor in contact 
with the stator; with the free body diagrams for each element provided in 
Figures 24(a) and (b). The equations of motion for the rotor and stator are 
developed according to the balance of forces in the free body diagrams. The 
me@ees) acting On the rotor are identical in nature to the forces present in 


Mae rigid Stator case. The balance of forces for the rotor is written: 


> - 
ai we) 7 
re hls. 


R(t) = (Ze) (27) 
The reaction force which acts at the rotor-stator interface though, is now a 
variable which depends on both the rotor's and stator's dynamic 
characteristics. Since the stator is non-rigid, the reaction force drives 
the stator out of line with the bearing center-rotor center axis. The 
resultant phase relationship between the rotor displacement and the rub 
location may be described by the angle A, which is a function of the rotor 
displacement, stator displacement and reaction force lines of action. 

Tne reaction force may thus be broxen down into its normal and tangential 


(friction) components with a periodic time dependence characterized by the 


synchronous speed, w, and the phase angle ): 


Gu Biel ae 90, 


R(t) = (N + juN)e 


or 
R(t) Na + 2) 172 QJ fut St 


where y = ee oa 
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FIGURE 2340). ROTOR-STATOR LOAD DIAGRAM (Top> 
FIGURE 2306), ROTOR-STATOR DISPLACEMENT VECTOR DIAGRAM 
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FIGURE 24. ROTOR AND STATOR FREE BODY DIAGRAMS 
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Additionally, the displacement of the rotor's center of gravity may be 


expressed: 
> > 


Ps 
Z, eee” . god cut 6) 


g 


8 is the phase angle which accounts for the lag between the rotor center 
meets forcing function. This forcing function now comprises both the 
eentrifugal force due to unbalance, Fos and the reaction force, R(t). The 
resulting expression for 6 is thus a function of three phase angles: 4, 
the phase angle due to rotor damping; w, the phase angle describing the 
menewor action Of the reaction force at the rub location; and i, the phase 
@rtterence between the rub location and the rotor center's displacement: 


9 = tan [uD _+ n(n)? sin(yeayI/E(K uM) PON Ca yy seca Chon 


Making these substitutions, the force balance equation for the rotor may 


now be written: 


| 2 | 
ee ee ce Sees es ged 2"), (085 
re Loe i at Cc 


Equation (28) may be solved for the displacement of the rotor center: 


172 juttr+p) 


um 6e9 OF" 8). 


2 
pal) 


Z 
L(K w M) + (jwD) J 
Breaking the rotor's impedance down to its magnitude and phase yields: 


+2) 1/2 QJ Cutt rty) 


ee e CE) 
2 2 se ae 
L(K - w M + (wD) ] 
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Dividing through by the term re and substituting non-dimensional 


variables, the rotor displacement may be finally written: 


18 2.1/2 jlwtra ; 
See en) ee Salem ey ae (29a) 
re (Ge & ean es 
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Whereas the centrifugal force due to unbalance drives the rotor's 
motlon, it is the synchronous rub force (i.e., the reaction force, R(t)) 
which acts as the input driving force for the stator's motion. The forces 
Pee respond to this driving force are the stator's inertial, elastic and 
damping responses. The force balance for the stator is derived from the 


free body diagram: 


eon Se oe Me (30) 
35 35 = M945 3 

. | = ie Sas tartan) 
By assuming a synchronous whirl response by the stator, Z = \z |e 


—> 
and solving for Z the stator's equation of motion is determined: 


2 > 
_ + 
[(K, W My) (jwD.)] 
Reducing the stator's receptance to a magnitude and phase, tT, and 


Substituting for the reaction force, the stator's equation of motion 


becomes: 


re ee ey eee oS nek 
= 2 2 Os 2 
Loe Ww M,) + (wD .) ] 


Maemwerc = Lan —---~5-- 


This equation may be expressed in non-dimensional terms by utilizing 


the following definitions: 


thes Agee 
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Now dividing through by the term ae anademaceoring Out the ratio (K 7K o)s 


the non-dimensional relation for the stator's displacement becomes: 


K : 7 
. (ary n(14 ye) 72 _5 YA 0) 
Z Ke 
a a era a coms endl ceramr cmerlgee are ae gad seen ae Nome iaeeei sor ier ese er eee C2 ) 
jJwt merce 22/2 
re C(1 B ys (27 8.) ] 


Two variables of interest in this analysis are the rotor's displacement 
and the reaction force at the rotor~-stator interface. To determine the 
rotor's displacement ae Equations (27) and (30) are utilized. The reaction 
force acting upon the stator is equal in magnitude but opposite in direction 
to that acting upon the rotor, so Equations (27) and (30) may be combined 


with the force term R(t) eliminated: 


> a > > 
pies h Zee? -=DoZae = Mz. + Mz (32) 
aC Seow co ses 
Figure 23(b) provides the vector triangle relating the stator displacement, 
rotor displacement and radial clearance. The kinematic condition which 


erases from the geometry presented here is given by: 


2, We eee, 
C S 


The implication here is twofold: 
(i) The magnitude of the rotor displacement relative to the stator 
displacement equals the radial clearance between the two elements. 
(ii) The point of contact between the two elements must be in line with 
the centers of the displaced rotor and stator. 
Applying this kinematic relationship and the previously derived condition, 
a4 god i wttd) 


> 
= As 
g C 


IeOLor 3; 


, results in the following equation of motion for the 


es r) Z 4 
2 > 3 > d soars o) jC wt+rA) 

“Pp —, entetiinmntitaall 
Be Z,* KD + Doz + (K+ Kidz. HE ee Sa 


imessyncechronous whirl solution is thus given by: 
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BOK K,) Ww (M., M,) J jlw(D, ae 
The time-independent, non-dimensional solution for the rotor displacement 


is obtained by dividing each side by the quantity eee and factoring out 


K on the right hand side. The resulting expression is: 


Z gel? + K el Ly 
RerQCMNEGEIEES§ "99 6G5°5° °° 9 500 ene nr ese eS (33a) 
re : K 5 ye 1, ? 12 
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This equation is further reduced to a single magnitude and phase in the 


following: 
K mK 
Z [ (8° cose 7 ==cosh)* : (ep°sing + =2sina)@]'/* 

c r . 
Rat oo Du us in a ra a K) (33b) 
ee @: Mo > Wh  _F 4 Bz 

VieGieem e+ =-( 1-8 [a a + [2c8 + 2cB ---] } 
A Ong Wy 6 
K 
(ep“sing + g>sinh) 
where v = ae ee re ee 
Cap aeen + z=cosh) 
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ime reaction force is more difficult to solve than the rotor displace- 
ment, but provides much meaningful information regarding the rub motion. A 
positive reaction force indicates full annular synchronous rub is in 
equilibrium. Conversely, a negative reaction force indicates a no-rub 
action, while an imaginary reaction force indicates full annular synchronous 


rub is not in equilibrium (i.e., not a solution to the rub equation). In 
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order to solve for the reaction force as a function of the rotor speed, the 
kinematic condition relating the rotor and stator motion must once again be 


agp lied: 


The rotor's displacement relative to the stator is derived from Equations 


(oa) and (31): 


K 
+ 2 4( 0-4) . (sF)e J" 
eres )e* | 5 __, 
n glut we 4 We Be 
c Je ey S 
where R= neiee oy + (278)°] 
Pe 2 
and Re = [(1 Bs ye (26 .8,) ] 


Applying the above kinematic condition, this equation may be written: 


(oS J Wat) 
ee a peeidc.c \ ie Ss JIA _ oA 
R ee R iw2 H We 
ié ie Ss 
melving for the phasor ev yields: 
/ ye 
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By defining 





-65- 


K 
Wes za es. eae 
r Samet rat Cia) 
K 


le ellie A eie Oo) + (SF yR | ’*cos( y-1)1 
S 


Care ee nila) LR 


then the relationship becomes: 
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Therefore, the phase angle may be defined in terms of known variables: 
Meri o= ) = 0) ; 

with the magnitude of Equation (32) being equal to unity. The square of 

the magnitude may also be equated to unity, yielding: 


K 


R (eB) /(E(RR_) + n(1eye)'/*cR.'/cos(y-9) + GIR, Peos(v-0) 1" 
D1 JS We Spy, 1/2 2 
+ inn)! *TR'sin(ye) + (ZF )R'/*sin(y-1)]}*) = 1 (3ua) 
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Equation (34a) is resolved into a quadratic equation in the variable n in 


the following expression: 


2 Be 1/2 K RAS ie : 
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Solving this equation for the normal force, n, yields the lengthy 


expression: 
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An alternative expression is derived by applying the following 


relationships: 
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It should be noted that the expression above reduces to the earlier 
derived result for rub against a rigid stator (Equation 11)) by employing 


the corresponding conditions: 


K = ©, therefore K /K_ = 0 
Ss rs 
NS 7 
- = 0 
B = 0 
R Sy 
= 


The impact of releasing restrictions on the stator's dynamic behavior 
is particularly complex due to the numerous variables involved. By 
comparing Equations (11) and (35) though, qualitative conclusions may be 
drawn. There are three dynamic ‘elements' significant to this system: the 
rotor, the stator and the combined rotor-stator. The character of the 
response of each element is dependent on the speed of excitation and the 
element's natural frequency. By considering the magnitude of a generic 
element's mechanical impedance, [(K - aie # Gye then it is apparent 
that at low excitation frequencies, the response is as a stiff spring, 
characterized by the value of K. At high frequencies, the response is aS a 
massive element; characterized by the value of M. The response iS a maximum 
in the vicinity of the natural frequency, V/K/M, (specifically at 


ene yomy!/2). 


W = Lu, These simple principles may be applied to help explain 
the dynamics of the rotor-stator system. 
Given the rotor's and stator's characteristic mass and stiffness; 


the composite element's natural frequency, Wh occurs between Onin and Ons? 





—-§8- 


+ K 
oer) 22.5 
“N 
r S 


ieee simple relationships thus exist between the rotor and stator: 


(i) the stator's natural frequency iS greater than the rotor's natural 
frequency 
(ii) the stator's natural frequency is less than the rotor's natural 
frequency 
(iii) the stator's natural frequency is of the order of the rotor's 


natural frequency 


Bach of these relationships may be satisfied by several conditions relating 
Eie elements’ mass and stiffness. For example, Wh > Wn is unconditionally 
true when K > K_andM_=M_or when K_ = K. andM <M. It is condi- 

Ss a — r S eS Ie 
meomally true wnen K > K andM »>M _ orwhenK <K_andM <M. It is 

Ss > Ss ie iS i ——= 5 ie 

evident that a large number of cases may be analyzed based on the relative 
masses and stiffnesses alone. This is exacerbated when variations in 
unbalance, friction and damping are investigated. 


iw 


Ns ls sufficiently greater than w 


Nr? then when rub first initiates, 
PiematalOr response resembles the infinitely stiff case; with the resulting 
reaction force depending directly on K.- In general, since it is the 
stator's stiffness which keeps it pressed against the precessing rotor at 
Speeds below Wes then continuous synchronous rub becomes a less stable mode 
of motion as Ke decreases. As rotor speed approaches Wy» the rotorssvator 
system becomes very compliant and continuous synchronous rub will cease 
(i.e., the radical term in Equation (36) becomes negative). As the rotor 


speed further increases and reaches w the stator alone is extremely 


Ns’ 
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compliant. Since a small push by the rotor thus brings a disproportionately 
large response by the stator, then continuous contact may rapidly diminish. 
(This is evident in Equation (36) since the value of RS tends toward zero as 
a” approaches 1.0 and the negative term within the radical, 

—(K /K.)°(R/Ro), becomes dominant). Throughout this speed regime, the 
rotor's motion (if unperturbed) is identical to that of the free rotor. At 
yet higher speeds, the stator's mass dominates the character of its 
response. Using limit analysis, the stator's behavior once again approaches 
the rigid response at very high speeds (8 >> Bas The criterion for 
Stability of synchronous rub at very high speeds is identical to that for an 


Merinictely stiff stator; i.e., 


then synchronous rub is stable for 8 >> Be 

At this speed though, the rotor must be perturbed to initiate the rub. 

ig Wiis and Wn are equal, then the stator's response resembles the 
rotor's response for any particular excitation frequency. Therefore, the 
Speed at which the combined rotor-stator is most compliant coincides with 
miewopeed at which the free-rotor's runout is greatest. This similarity in 
Pesponse tends to nullify the stator's non-rigid character if (K/ Ko) and 
(M/M.) are greater than 1.0. As Ke and My are increased, the response 
rapidly approaches the rigid stator case. When the stator's dynamic 
Characteristics equal the rotor's, then the speed regime in which continuous 
Syoenronous rub is a solution to rotor motion is identical to the rigid 
Stator case and the reaction force is exactly half that which would result 


if the stator were rigid. As the stator becomes less stiff and less 


massive, synchronous rub becomes less stable. In the limit, for low values 
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of (Ko/K) the system is unable to sustain a continuous rub with the 
precessing rotor. 

If the stator's natural frequency is measurably less than the rotor's, 
then the speed regime in which the stator is most compliant occurs prior to 
Mreminitiation of rub. When the rotor first gains contact, the stator's 
impedance and that of the combined rotor-stator are increasing with the 
Bauare Of the speed. The response of the system rapidly approaches that of 


the rigid stator case as Me increases and thus either w decreases or, for 


Ns 


a constant Wag? K, increases. If in fact a low value of Wr is due to low 
values of both Ke and Mo: then On may be approximated by Or and the rub 
response is as described previously for Oy = Whe 

inewimpact Ol the Stater’sS damping characteristics are straigntforward. 
Independent of the natural frequencies, stiffnesses or masses of the rotor 
and stator; any increase in the stator's damping coefficient tends to 
diminish the equilibrium of synchronous rub. Assuming light damping for 
both rotor and stator, variations in stator damping has little effect on the 
magnitude of the reaction force during rub. 

In summary, if the stator's natural frequency iS greater than the 
mecor S Natural frequency, then synchronous rub is not in Sour 
throughout a speed regime which encompasses the stator's and the combined 
rotor-stator's natural frequencies Cw and Wy respectively). This speed 


regime widens for a given value of w aS the stator stiffness, Kos and 


Ns 
Stator mass, Moo are decreased. Similarly, the reaction force experienced 
during continuous rub for a given value of Ws decreases as K and M. are 
decreased. As the stator's natural frequency is reduced towards the rotor's 


natural frequency, this particular speed regime where synchronous rub is not 


Viable diminishes accordingly. When Wig = Wr? the stability of synchronous 
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rub resembles that of the rotor against a rigid stator; being more stable if 
(K/K.) (and thus (M/M)) is less than unity and vice versa. As may be 
expected, the reaction force in this case increases as Ke and ie iMmCcredae. 
As the values of Us, and Wis are reduced below Wr the synchronous rub system 
response becomes as stable as if the rotor were rubbing against a rigid 
Stator. This is true since the stator's impedance rapidly increases with 
the square of the excitation frequency. Regarding stator damping, any 
increase in the damping coefficienct, Do» reduces the speed regime in which 
synchronous rub may occur. In general, the stator's finite mass, stiffness 
and damping characteristics serve to abate the negative characteristics of 
rotor rub. In no case do these characteristics intensify the reaction 
forces between the two elements (when compared with rub against a rigid 
stator), but rather they tend to reduce the reaction forces, often 
sufficient enough to create a dynamic relationship less favorable for full 
annular synchronous rub (and too, reverse whirl). 

Appendix A presents a numerical investigation of the impact of varying 
the parameters discussed in this section. Through the efforts of Professor 
Udo Fischer of the University of Magdeburg, Magdeburg, East Germany, a 
Fortran program has been developed which clearly demonstrates the time 
@epenaent response of the rotor and stator in accordance with the equations 
of motion developed here. The stability of the synchronous rub precession 
in those speed regimes where this mode of rotor motion is in equilibrium is 


clearly demonstrated. 





ao 


6. CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 


Mechanical vibrations in rotating machinery result from any mass 
unbalance present in the rotating element. Should rotor rub occur as a 
result of these vibrations, the integrity of the system may be significantly 


impaired. Two direct solutions to this problem are immediately apparent: 


(i) Reduce the mass unbalance which drives the vibrations in the 
eyocem, or 
(ii) Increase the radial clearance between the rotating element and the 


Stationary element. 


Both of these solutions are costly, either through the precision required in 
the manufacturing phase or efficiency lost throughout the operational life 
of the system. 

The dynamics of rotor rub are quite complex, but in general a system 
may be resolved to one of four behavioral types dependent on the system's 


parameters: 


(i) Synchronous precession without rub 
Sr Partial rotor rub 
is) Full annular rotor rub 


(iv) Reverse whirl 


Most rotating machinery are designed for the no-rub condition. Analyti- 


cally, rub may be anticipated if the ratio of the mass eccentricity~to-the- 
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radial clearance exceeds twice the damping ratio for the system. Should rub 
eecurae tie equllibrium of full annular rub 1s then dependent on both the 
rotor's and stator's dynamic characteristics, the system's unbalance, the 
radial clearance and friction experienced between the two contacting 
surfaces, and the rotor's speed. 

Initial analysis was conducted with the assumption that the system's 
Stator is infinitely rigid. Dependent upon the system's damping ratio, 
displacement-to-clearance ratio and frictional coefficient; the rotor's 
response as a function of speed can then be characterized according to one 


of the following: 


Semeinited syionronous rub: synchronous rub is in equilibrium up to a 
"critical speed" (unique to each system) above which equilibrium is never 
regained. 

Sey eoynenronous rub = limited discontinuity: Similar to the limited 
Eupecase, equilibrium is lost at some critical speed, but is then regained 
at a higher speed. 

(iii) Continuous synchronous rub: As the description implies, this type 


System 1S able to find an equilibrium position for all rotor speeds. 


Speed regimes where full annular rub is either not in equilibrium or 
not stable may fall into a partial rub (if this mode is in stable equili- 
brium at this rotor speed) or jump down to a synchronous whirl without rub. 
An alternative mode of rotor motion which may be stable under these 
conditions is a reverse whirl driven by the counter-torque produced by the 
frictional forces acting on the rotor. Some disturbance is necessary to 


perturb the system into the reverse whirl, but such may occur through a 
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Meeotene sacl ial rub motton following the deterioration of full annular 
IUD. 

Expanding the analysis to investigate the impact of the stator being 
non-rigid results in extremely detailed equations describing the motion of 
the system and the forces at work between the rotating and stationary 
=eenenc. in contrast to the rigid=stator system, the stator now joins with 
the rotor in responding to the unbalance driving the vibrations. In 
general, equilibrium of the full annular rub diminishes in the vicinity of 
the combined rotor-stator natural frequency (which previously was infinite). 
Additionally, the amplitude of the motion varies greatly now that the stator 
is not absolutely constraining. If the stator natural frequency is greater 
than the rotor natural frequency, then the recurrence of rub at speeds above 
pos requires some perturbation to initiate contact. In general, the finite 
dynamic characteristics assigned to the stator which make it non-rigid 
diminish the reaction forces experienced by the rotor and stator (as 
compared to the rigid stator case) and in fact may create a dynamic 
merarrousnip less favorable for full annular synchronous rub. 

Much work remains to be accomplished in the study of rotor rub. 
Experimental investigation of the results presented here would prove 
particularly meaningful in validating the results and justifying the 
assumptions made. An analytical investigation of the equilibrium conditions 
for partial rotor rub would be very complementary to this study on the full 
annular rub problem. Of particular significance would be work which 
contributes to the analytical determination of the stability of full annular 
rub, partial rub and reverse whirl for a rotational system as a function of 


the excitation frequency. 
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APPENDIX A 


It is the intent in this Appendix to investigate the dynamics of 
several exemplary rotor-stator systems in accordance with the analysis thus 
mareconducted for rotor-rub. Unlimited combinations of rotor and stator 
dynamic characteristics are possible, but here a limited number are studied 
to demonstrate the impact on system behavior of varying key parameters. In 
particular, variations in rotor and stator natural frequencies, stiffnesses, 
masses and damping characteristics are investigated. Similarly, the effects 
of varying the magnitude of the unbalance and friction are studied. The 
indicators for system behavior deemed here to convey the most information 


are: 


(i) Rotor Amplitude vs. Speed 
(ii) Normal Force vs. Speed 


(iii) Rotor and Stator Displacement History 


Synchronous rub is presumed in the analysis, and is indicated by stable 
rotor displacement solutions (aside from the fundamental free-rotor 
solution) in conjunction with a real, positive reaction (or normal) force at 
the rotor-stator interface. 

The plot of rotor amplitude (Zz /r 4) vs. frequency ratio (w/w? is 
presented for the two solutions to Equation (29) as well as the free-rotor 
Solution to Equation (4). The free-rotor solution (depicted by the broken- 


line) is the only stable solution below the rub speed, 8 and is again 


rub’ 


Stable at speeds greater than 8' (described in section 3). Both equilibrium 


solutions to Equation (29) are presented (depicted by the solid lines 








70> 


iapelicd leand 2). but 1t 1S only curve 1, that solution initiating at the 
rub speed, which provides stable equilibrium. This is somewhat intuitive by 
analyzing the geometry of the rotor and stator equilibrium conditions (such 
as Figure 23) corresponding to a particular speed of interest. Where the 
solution is discontinuous, the recurrence of this stability is not assured. 

The normal force is also plotted for the stable solution as a function 
of rotor speed. Regions where the normal force is negative indicates rub is 
feelin action. IU iS immediately apparent that no solution exists for 
eyaenronous rub in the vicinity of the combined rotor-stator and the stator 
natural frequencies. Stability of equilibrium is determinable using a 
Simulation program provided by Professor Udo Fischer which traces both the 
rotor's and stator's displacement as a function of time. Several examples 
emepotn Stable and unstable equilibrium are provided throughout this 
Bppendix. A Stable equilibrium is indicated by the constant orbit of both 
the rotor and stator. This solution may be artificially perturbed (through 
initial conditions to the program), but still returns to the equilibrium 
condition. Unstable equilibrium is evidenced by varying amplitude of both 
rotor and stator displacement, typically ending in the rotor undergoing a 
Seyucnronous whirl without rub while the Stator returns to the rest position. 

Two primary groupings are analyzed: Ws 2 Oe and Ws < Wir? A total 
of 18 cases are presented with their key parameters outlined in Table Al. 
Rotor characteristics are indicated by subscript 1 while stator characteris- 
ties are indicated by subscript 2. One baseline case for each of the 
primary groupings is analyzed, with comments regarding each subsequent case 
within the grouping being referenced to the baSeline case (unless specifi- 
cally stated otherwise). 

These cases, as described in the following, clearly demonstrate the 


impact of varying the system's parameters as previously discussed: 





a Oe 


Case 1: This case provides the baseline for the condition igs > Wn As an 
aiprosor limited Synchronous rub, Only a limited speed regimes exists 
wherein synchronous rub is in equilibrium. 

Case 2: This case demonstrates the impact of an increase in the mass 
unbalance. AS an example of synchronous rub-limited mass unbalance. As an 
eramole Of Synchronous rub-limited discontinuity, synchronous rub is in 
equilibrium for all speeds greater than the rub speed, except for a limited 
region demarked by the combined rotor-stator and the stator natural frequen- 
cies. Additionally, the amplitude of displacement is increased above the 
baseline while the reaction force increases severely. 

Case 3: A reduction in friction increases the speed regime wherein 
synchronous rub is in equilibrium. Concurrent with this is a large increase 
mathe reaction force. 

Case 4: Similar to Case 3, a reduction in stator damping leads to estab- 
lishing the equilibrium of synchronous rub over all speeds outside of the 
combined rotor-stator and stator natural frequencies. As well, there is an 
increase in the reaction force. 

Case 5: Reducing the rotor's damping results in a large increase in the 
amplitude of displacement (above the baseline) and a slight increase in the 
Speed region for synchronous rub. No recurrence of rub occurs though above 
the stator natural frequency. 

Case 6: An increase in excitation and damping on the rotor and stator 
demonstrates how the damping reduces the displacement and reaction force at 
low speeds (in comparison to Case 2), but has little impact at high speeds 
(w > ug): 

Case 7: Reducing stator rigidity while maintaining the same stator natural 
frequency results in a reduced region of synchronous rub and a reduced 


reaction force. 
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Case 8: Converse to Case 7, increasing stator rigidity increases the region 
of synchronous rub and results in increased reaction forces. 

Case 9: Increasing stator stiffness and so the stator natural frequency 
extends the speed regime for synchronous rub while cauSing an increase in 
the reaction force and a reduction in displacement. 

Case 10: Increasing stator stiffness and the excitation demonstrates the 
large displacements and severe reaction forces that may develop. Once 
again, synchronous rub is re-established above the stator natural frequency. 
Case 11: This case represents the baseline for the condition Wis < One A 
limited synchronous rub occurs with a displacement less than the radial 
clearance. 

Meeemlc: increasing excitation results in a continuous synchronous rub with 
Significant reaction forces developing at high speeds. 

Case 13: Reducing friction also causes a continuous synchronous rub to 
occur similar to the previous case. 

Case 14: Reducing the stator's damping has only a slight impact on the 
reaction force and equilibrium range for synchronous rub at speeds much 
greater than the stator natural frequency. 

Case 15: Reducing rotor damping results in a continuous synchronous rub 
Similar to Case 12. 

Case 16: Increasing stator rigidity, while maintaining constant natural 
frequency reduces the speed regime of synchronous rub in this case while 
causing the amplitude to approach the radial clearance. 

Case 17: Converse to the previous case reducing stator rigidity for a given 
Wr» causes a reduction in amplitude with a slight increase in the speed 


regime for synchronous rub. 
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Case 18: Increasing stator mass (and thus reducing uy ) has little effect 
on the speed regime or reaction force of the synchronous rub, but does drive 
the displacement toward the radial clearance (i.e., the rigid stator 
eondition). 

The stability of the solutions represented by displacement curves 1 and 
ees analyzed for the three rub behavior types: limited synchronous rub 
(Cases 1 and 11), synchronous rub-limited discontinuity (Case 2) and 
continuous synchronous rub (Case 12). Within these examples, several speeds 
are investigated in order to study the stability of each speed regime 
wherein synchronous rub is in equilibrium and too the stability of the two 
solutions at those speeds. These examples demonstrate the inherent 
instability of curve 2 solutions as well as the variability of stability 
associated with curve i solutions. In particular, curve 1 appears stable 
morepoth the limited rub and continuous rub conditions, but there is a speed 
range within the synchronous rub-limited discontinuity equilibrium solutions 


wherein Stability is lost. 
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TABLE Al 

Case No. K2/K1 M2/M1 D1 D2 
1 4.0 ee O.2 OZ 
2 4.0 en 8, Oz2 Ore 
3 4.0 ire 0 Gre Ome 
y 4.0 lene OR On.0 
5 4.0 (ale ‘Oa ie Gee 
6 4.0 e0 0.4 0.4 
fe ae Cn25 Ore Oe 
8 8.0 20 OU Oe2 
9 imomae) 1 0.2 One 
10 16.0 1.0 2 Oa 
11 Hale. 4.0 Ox2 0.2 
2 lee 4.0 Ore OLE 
13 152.0 4.0 m2 OR e 
14 10 4.0 Or 2 CAG 
ie ale, 4.0 O.0 One 
16 BAN) S20 One Bie 
int O225 ee) One OF:2 
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